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Abstract 

This article deals with different generalizations of the discrete stability property. Three 
possible definitions of discrete stability are introduced, followed by a study of some par¬ 
ticular cases of discrete stable distributions and their properties. 


1 Introduction 


Stability in probability theory refers to a property of probability distributions when a sum 
of normalized, independent and identically distributed (i.i.d.) random variables has the 
same distribution (up to scale and shift) no matter how many summands we consider. 
Random variables with this property are called stable and they form a wide class of 
probability distributions. Except for one particular case, the Gaussian distribution, all 
stable distributions are heavy tailed. The classical stability refers to stability under sum¬ 
mation but the concept can be extended onto other systems as well. Stability under 
maxima (or max-stability) leads to heavy tailed distributions called generalized extreme 
value distributions; stability under random summation where the number of summands 
is a random variable leads to heavy tailed j^-stable distributions. Stability of discrete 
systems is a topic that has not been studied as extensively as others but here also the 
discrete stable distributions exhi bit heavy ta ils. 

Introduced by Paul Levy in (iLevvl 11925!) , stable distributions are a generalization of 
Gaussian distribu tion in sever al wa ys. The theory of stable distributions was developed 
in monographs by Le^ j[l9^ ) and Khintchin3 ( 19^ h 8-nd further extended in the work 
by iGnedenko and Kolmogorw ( 19491) anT Feller ( 19701) . There exist few equivalent defi¬ 
nitions of stable distributions. Paul Levy defined stable distributions by specifying their 
characteristic function. For that he used the Levy-Khintchine representation of infinitely 
divisible distributions. Second definition is connected to the “stability” property - a sum 
of stable random variables is again a stable random variable, a well known property of 
Gaussian random variables. Third is the generalized central limit theorem - stable distri¬ 
butions appear as a limit of sums of independent and identically distributed random vari¬ 
ables without the standard assumption of the central limit th eorem about finite variance. 
This r esult generalizes the central limit theorem and is due to IGnedenko and Kolmogorov! 
( 19491) . Gaussian distribution is a special (limit) case of stable distributions, the only 
stable law with finite variance. Recent and extensive overview of the theory of stable 
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random variables can be found i n Zolotarev ( 19861) . Uchaikin and Zolotarev ( 1999f) and 
Samorodnitskv and TaaoJ ( 1994h . 

In many practical applications continuous distributions are often preferred over dis¬ 
crete distributions because they offer more flexibility. There are however cases of practical 
applications where one ne ed to describ e heavy tails in disc rete data. C itations of scientific 
papers (first observed by Pricel ( IQbS P. word frequency ( Zipi ( 1949ll ') and population of 
cities are all well known examples of discr ete d a ta wi th power tails. A simple discrete 
power law distribution was introduced by H (I 1949 I 1 and relied on the zeta function 
(therefore called Zipf or zeta distribution). 

Another possibility is to consider discrete variants of stable and ^-stable distributions. 
The notion of discrete stability for lattice random variables on non-negative integers was 
introduced in ISteutel and van Harn ( 19791) . They introduced so called binomial thinning 
operator 0 for normalization of discrete random variables. That means that instead of 
standard normalization aX by a constant a S (0,1), they consider aQX — where 

ti are i.i.d. random variables with Bernoulli distribution with parameter a. As opposed 
to the standard normalization, this thinning operation conserves the integral property 
of a discrete random variable X. Together with a study of discrete self-decomposability 
they obtained the form of generating function of such discrete stable distributions. By 
considering only non-negative discrete random variables, they obtained a discrete version 



probabilities as well as their asymptotic behaviour. They showed that the discrete stable 
distribution belongs to the domain of normal attraction of stable distribution totally 
skewed to the right with index of stability smaller than one. The non-existence of a closed 
form formula of the probability mass function and non-existence of moments implies that 
the classical parameter estim ation procedures such as maxi mum likeliho o d and method of 
moments cannot be applied. Marcheselli et al. ( 2008 ) and Dorav et al. ( 20091) suggested 
some methods of parameter estimation of the discrete stable family based on the empirical 
characteristic function or on the empirical probability generating function. 

Discrete stable distributions in limit sense on the set of all integers were introduced 
in Klebanov and Slamov^ ( 201,11) . Two new classes of discrete distributions were intro¬ 


duced , generalizing the definition of discrete stable distribution of ISteutel and van Harn 
( I 979 I) on random variables on the set of all integers. It was shown that the newly in¬ 
troduced symmetric discrete stable distribution can be considered a discrete analogy of 
symmetric a-stable distribution with index of stability a G (0,2], whereas the introduced 
discrete stable distribution for random variables on 1 can be viewed as a discrete anal- 
og v of g-stable distribution with index of stability a G (0,1)U{2} and with skewness 
/3. ISlamova and KlebanovI (120121) gave two distributional identities for the symmetric dis¬ 
crete stable and discrete stable random variables, allowing for simple random generator. 
Possible estimation procedures for the class of discrete stable laws were also considered. 

The aim of this paper is to study different generalizations of the strict stability prop¬ 
erty with a particular focus on discrete distributions with some form of stability prop- 
erty. The starting p o int of the article are discrete stable distributions introduced in 
Steutel and van HarnI ( 1979|) . Their definition of discrete stability is a simple general¬ 


ization of the classical stability property where they consider only one type of thinning 
operator. The classical stability property can be formulated in several equivalent ways 
and our aim is to study generalizations of these equivalent definitions for the discrete 
case. We propose three definitions of discrete stability for random variables on non¬ 
negative integers. The main focus is on the first definition that generalizes the definition 
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of Steutel and van Harn ( 19791) by allowing the thinning operator to be an arbitrary dis¬ 
tribution satisfying certain condition. We introduce also the symmetric and asymmetric 
variant of disc rete stable distribution. The d efinition of discrete stability on all integers, 
similarly as in iKlebanov and Slamovai ( 2013 1. is possible only in the limit sense. 

In Chapter 2 three possible definitions of discrete stability for non-negative integer¬ 
valued random variables are given. These definitions consider different approaches to 
introducing discrete stability, each of them being a discrete version of a different definition 
of stability in the usual sen se. The first definition generalizes the approach taken by 
ISteiitel and van HarnI ( 1979ll and considers a general thinning operator to normalize the 
sum of discrete random variables. The second definition takes the opposite path and 
uses a general so called portlying operator to normalize discrete random variables. The 
last approach combines the two definitions and as it turns out includes the previous two 
definitions. Examples of the thinning and portlying operators for which a positive discrete 
stable random variable exists are provided. Chapters 3 to 6 are dedicated to the study of 
analytical properties of discrete stable distributions in the first sense. The study is focused 
mainly on the class of distributions connected to modified geometric thinning operator 
and we give results on characterizations, probabilities, moments, limiting distributions 
and asymptotic behaviour for positive and symmetric discrete stable random variables. 
Section 6 gives also some results on properties of positive discrete stable random variables 
with Chebyshev thinning operator. 


2 On definitions of discrete stability 


ISlamova and Klebanov! ( 014bh introduced a possible approach to obtain discrete analogies 
of stable distributions. By approximation of the characteristic function of stable distribu¬ 
tion or of its Levy measure three discrete distributions were obtained. These distributions 
are discrete approximations of the stable distributions and it is not clear what properties 
they share with the stable distributions - by the construction it is obvious they have the 
same tail behaviour, but it is not clear whether they share other properties as the stabil¬ 
ity property, self-similarity, infinite divisibility and others. In this Section we define three 
new classes of discrete probability distributions by generalizing the stability property for 
discrete random variables. 

The strict stability property of continuous random variables can be defined in several 
ways. We say that a random variable X is strictly stable if one of the following holds 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


n 

cX = aXi + bX 2 , 


where Xi, X 2 , ■ ■ ■, Xn are independent copies of X and a„, a,b and c are positive 
constants. If we want to define a discrete analogy of stability we have to reconsider the 
normalization by the constants a„. An, and a, b and c, as the normalized random variables 
are not necessarily integer-valued. We may consider the following modification. Consider 
for example the first definition and let us assume that X is non-negative integer-valued 
random variable. We may write 


X = l-|-1-1-- hi, and pX = p + p-\ - hp, 

X times X times 
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where we normalize X by a constant p G (0,1). Instead we can consider a thinning operator 
p Q X, where 

J 30 ^ = £ 1+£2 + ‘‘‘ + ) 

'-V-' 

X times 

where are i.i.d. Bernoulli random variables with Ee^ = p, i.e. 

__ J 1, with probability p, 

* \ Oj with probability 1 — p. 

In the following Sections we introduce three different definitions of discrete stability 
generalizing the definitions of strict stability (1^ ~ (EHI for the case of non-negative 
integer-valued random variables. 

2.1 On first definition of discrete stable distributions 

In this Section we give a definition of discrete stability that generalizes the first definition 
of strict stability (12.11) for discrete random variables. The multiplication by a constant a„ 
can be understood as a normalization of the sum Xi, or normalization of the individual 
summands Xi. In the case of discrete random variables one can not use this normaliza¬ 
tion as it violates the integral property of the summands. We need to find a different 
normalization that maintains the integral property. One possibility is to use the binomial 
thinning operator 


X 

X{a) = aQ X = ^ €i, where F{ei = 1) = 1 — F{ei = 0) = a, 

i=l 


instead of QnXi. This normalization was used in Steutel and van HarnI ( 1979ll to define 
discrete stability on Nq. One can generalize this definition of discrete stability by consid¬ 
ering a general normalization, or “thinning” operator. 


Definition 2.1. Let X, Xi,X 2 , ■ ■ ■, X^ ... denote a sequence of independent and iden¬ 
tically distributed (i.i.d.) non-negative integer-valued random variables. Assume that for 
every n G N there exists a constant Pn G (0,1) such that 


n Xi 

(2.4) X = '^Xi{pn), where Xi{pn) = Pn & X^ = '^£j'\pn), 

i=i j=i 

and £^j\pn) are i.i.d. non-negative integer-valued random variables. Then we say that X 
is positive discrete stable random variable in the first sense. 

This definition is rather general as it offers a flexibility on the choice of the “thinning” 
distribution of random variables e. This flexibility is however limited as a positive discrete 
stable random variable exists only for some choice of the thinning distribution. A question 
is therefore how to describe the family of thinning distributions for which a positive 
discrete stable random variables exists. 

Let us denote the probability generating functions of the random variables X and £(pn) 
by 'P{z) = E[z^] and Qp„iz) = E[z®(P")] respectively. There is an equivalent definition of 
positive discrete stability in terms of those probability generating functions. 

Proposition 2.2. A random variable X is positive discrete stable if and only if for all 
n G N there exists a constant Pn G (0,1) such that 


(2.5) 


p(z) = p-(Qp„(z)). 
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Proof. It follows from the definition (12.41) that X is positive discrete stable if and only if 

Viz) = [V^iz)]\ 

The probability generating function of X can be computed in the following way. 


r^{z) = E = k)E 


Y,nx = k) (e 




y<^l{Pn) 


-\ \ k 




= nQpAz))- 

Hence X is positive discrete stable if and only if its probability generating function satisfy 
the relation 

p(z) = P"(Qp„(z)). 

□ 


Remark 2.3. It follows from the definition that a positive discrete stable random variable 
X is infinitely divisible: for every n G N there exist random variables Yi,Y 2 ,... ,Yn such 
that 

X = Yi+Y2 + --- + Y„. 

This obviously holds for Yi = Xi(pn). 

Further denote by Q a semigroup generated by the family of probability generating 
functions {Q{z) = (z),n G N} with operation of superposition o. It can be shown that 

a superposition of two probability generating functions is again a probability generating 
function. 

Lemma 2.4. If Qi(z) and Q 2 iz) are two probability generating functions of two random 
variables with values in No, then their superposition 

QioQ2{z) := Qi{Q2{z)) 

is also a probability generating function of some random variable with values in Nq. 


Proof. Let be a random variable with values in Nq with probability generating function 
Qi and Xi,X 2 ,... i.i.d. random variables with values in No with probability generating 
function Q 2 . Define a new random variable S by 

N 

s = Y.x,. 

i=l 

Then S' is a random variable with values in No. Its probability generating function can 
be computed using the fundamental formula of conditional expectation as follows 


QsW=E[z''] =E 




= ^P(lV = n)E 


n—0 




= ^ F{N = n)E ^ P(fV = n) [Ez 




n—0 


n—0 


= J2nN = n) [Q 2 {z)r = Q^{Q 2 {z)). 


n—0 


So the superposition Qi o Q 2 {z) is a probability generating function of random variable 
S with values in No. □ 
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Now we show that the semigroup Q must be commutative. 

Theorem 2.5. Let X be a positive diserete stable random variable. Then the semigroup 
Q must be commutative. 

Proof. Let us denote G{z) = log'P{z). Then (12.511 is equivalent to 
(2.6) G(z) =nG(Qp„(z)), n G N. 

Let G(z) be a solution of (12.61) . Then for all n G N it must hold 

Qpniz)=G-^ 

It follows from here that for all ni, n 2 G N 

= G-1 (——G{z) 

\ni n2 
“ Qpna {Qpni (^)) ) 

which means that Q is commutative. □ 


Similarly as for the classical stable distribution, we can show that the constants 
have to take form for some 7 > 0 . 

Theorem 2.6. Let X be a positive discrete stable random variable in the first sense. 
Then there exists 7 > 0 such that pn in takes form 


Pn = n 


Proof. The proof follows ( Uchaikin and Zolotarev , 19991 §2.4) where a similar statement 
for stable distributions is proved. From the definition it follows that for every n > 2 we 
have X = ^i{Pn) where Xi,X 2 ,... are independent copies of X. Then 


X = P2 & Xi + P2 & X 2 , 


therefore also 


X = P2 Q {p2 © Xi + P2 0 ^ 2 ) + P2 © {P2 0 + P2 © Xif). 

But the operation 0 is associative: p Q {p Q X) = p'^ Q X. Let us denote Y = p 2 Q Xi + 
P 2 © X 2 . Then (using result from proof of Proposition 12.211 

'Pp.eviz) = Vy{Qp,{z)) = V\QpAQpM)) = 
because Q is commutative. Therefore 

X =pIqXi+pIqX2+pIqX3+pIoX4 


and similarly for every n = 2 ^ 

(2.7) X =p^QX4+p^QX2 + ---+p^QXn. 
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On the other hand, we have 


(2.8) X — Pn Q Xi + PnG) X 2 + • • • + Pn 0 Xn- 

Comparing (12.711 with (12.8L with n = 2^, we have Pn = P 2 - Hence 

logp„ = fclogp 2 = 7 ^ 1 ogP 2 = 
log 2 

So we obtain that 

Pn = n~^l^'^, 72 =-log 2 /logp 2 > 0 , n = 2 '‘,/c = 1 , 2 ,.... 

In a similar way, starting with sums with 3 terms X = p^Q Xi + pa 0 X 2 + pa 0 X 3 , we 
get 

p„ = n“^/'^^ 7 a =-logS/logpa > 0 , n = A: = 1 , 2 ,.... 

And in general case, 

Pn = 7 m =-log to/ log Pm > 0, n = m^,k = l,2,.... 

But for TO = 4 we obtain both 74 = — log 4/log p 4 and logp 4 = — 1/72 log 4. Hence 74 = 72 - 
By induction we conclude that 7 ^ = 7 for all to and therefore 

Pn = , for all n > 2. 


□ 

The question is how to extend the definition of discrete stability to contain not only 
random variables on No, but also on the whole integers Z. It is obvious that the sum in 
definition of X does not make sense for random variables that can achieve negative values. 
One possibility is to take the positive and negative part of X separately and consider again 
the same thinning operator. We can, however, obtain a wider class of distributions if we 
assume a different thinning operator than in Definition 12.11 

Definition 2.7. Let A, Ai, A 2 ,..., A^,... denote a sequence of independent and iden¬ 
tically distributed (i.i.d.) integer-valued random variables. Assume that for every n S N 
there exists a constant p„ G (0,1) such that 

n X+ X- 

(2.9) A = lim ^ Ai(p„), where W(p„) = ^ 4*^ (p„) - ^ (p„), 

i=l j=l j=l 

s^j\pn), €^j \pn) are i.i.d. integer-valued random variables, and A+ and X~ are the pos¬ 
itive and negative part of A, respectively (i.e. A+ = A if A >0 and 0 otherwise, 
X~ = —A if A < 0 and 0 otherwise). Then we say that A is discrete stable random 
variable in the limit sense. 

The main difference is that we do not assume the random variables e, e to be non¬ 
negative. The definition of discrete stability is only in the limit sense, not the algebraic 
one where we have equivalence in distribution in (EH) instead of the limit. 

Let us denote again the probability generating function of the random variables A and 
e{pn)) (and also e{pn)) by V{z) = E[z^] and TZp^{z) = respectively. 

We denote by Vi the generating function of the sequence {at = IP(A = fc), A: = 1,2,...} 
and by 7^2 the generating function of the sequence {bk = IP(A = k),k = —1, —2,... }. We 
denote Vq = P(A = 0). It is obvious that the generating function of A+ is Vq + Viiz), 
and the generating function of A“ is 7 ^ 2 ( 2 ). There is an equivalent definition of discrete 
stability in the limit sense in terms of those generating functions. 
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Proposition 2.8. A random variable X is discrete stable in the limit sense if and only 
if for all n G N there exists a constant pn G (0,1) such that 

(2.10) P(z)= lim [Po + Pl(7^p„(^))+P2(7^p„(l/z))]”. 

n—¥(X) 

Proof. It follows from the definition (12.91) that X is discrete stable if and only if 

V{z)= lim [V^izr. 


The probability generating function of X can be computed in the following way. 


= E P(X = fc)E 


k— — c 


\x = k 


Z^J 


oo ^ —1 

Y^{X = k)(E H ^{X = k)(E ) 


— k 


k—0 k— — oo 

= Po + ViiTZp^iz)) + V2(np^ (1/^)). 

Hence X is discrete stable if and only if its probability generating function satisfies the 
relation 

V{z)= lim [Po + Pl(7^p„(z)) + P2(7^p„(l/z))^. 


□ 


It is important to note that we do not define discrete stability property in the algebraic 
sense as we defined it for the non-negative integer-valued random variables. This also leads 
to the fact that we have no condition on the thinning operator TZ similar to Theorem 12.51 
In the following Subsections we introduce some examples of commutative semigroups 
Q leading to different positive discrete stable random variables. We will also give corre¬ 
sponding examples of discrete stable distributions in the limit sense. The proofs of the 
results will be provided in Chapter 5. 


Binomial thinning operator. Assume that the probability generating function Q 
is that of Bernoulli distribution with parameter^ G (0,1), i.e. we have Q{z) = pz+{l—p). 
It is easy to verify that the semigroup Q generated by probability generating functions of 
this form is commutative, as 


Qpi {Qp2iz)) =PlP2Z+ {1-P1P2). 

This operator was used in Steutel and van HarnI dlQTOll to define discrete stable distri¬ 
bution on No and it was showed there that it leads to a distribution with probability 
generating function given by 


(2.11) P(z) = exp{-A(l-z)^}, 7e(0,l], A>0. 


To obtain a generalization of this distribution on Z we can consider two-sided binomial 
thinning operator defined as TZ{z) = (1 — p) -I- pqz + p(l — q)z~^, where q G [0,1]. This 
thinning operator leads to a distribution on Z with probability generating function given 
by 


P(z) = exp 




IN'" 

1-qz- (l-q)-] -A 



1 

l-q--(l-q)z 


with A > 0,7 G (0,l],/3 G [—l,l],g G [0,1]. We can see that for /3 = 1 and q = 1 the 
distribution reduces to positive discrete stable (12.111) . 
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Thinning operator of geometric type. A generalization of the previous example 
can be obtained if we consider Q to be the probability generating function of modified 
geometric distribution with parameters p S (Oj 1) and /c G [0,1). Consider a function 


(2.12) Q(z) 


/ {l-p) + {p- k)z'^ \ ™ 
\ (1 — pn) — k{ 1 — p)z™ J 


{0 < K < 1, 0 < p < 1, m = 1} 
or 

{0 < p < K < 1, m G N, m > 1}. 


Lemma 2.9. The function Q(z) is a probability generating function. 

Proof. To verify that Q(z) = 9™-^" i® a probability generating function we have 

to show that the generating sequence {q„, n = 0, Ij • ■ ■} is a probability mass function, 
i.e. Y^^qn = 1, and 0 < (?„ < 1. We see that Yn^n = Q(l) = 1- We expand Q into 
a power series to obtain the generating series {qn,n = 0,1,...}. We will treat the case of 
m = 1 and m > 1 separately. 


Let first m = 1. Then we obtain Q(z) = Y’^=o with 


90 


1 -P 

1 — PK ’ 


qn=p ^ 


(1—p)” ^(1 — 

(1 — pk)"+i 


n > 1. 


We can easily verify that for 0 < k < 1 and 0<p<l, {g^jisa probability mass function 
and thus Q is a probability generating function. 


Let TO > 1. We obtain 


c(0 = H 


QnZ 


n—0 


where the coefficients q„ are given as 


<ln=Y. 

i=o 


1 \ l/m+n—7 / \ 7 /, / , ■ 1 \ /I / 

l-p\ ■' (p - K,y „_/l/TO + n-J-l\/l/ 


1 -p 


1 — PK 

This can be reduced to 

1 — p \ /\ j'jfYi n — 1 


qn = K 


1 — pK 


n 


n-j 


2 F 1 I {—1/to, —n}, 1 — 1/m — n. 


n G Nq. 


{k-p){1-Pk) 

(1 — pYk 


It follows from the properties of the hypergeometric 2 F 1 function that 0 < < 1 if and 

only if 

(/c-p)(l-p/c) 

- (1-p)2k - ■ 

This is fulfilled if and only if0<p< k < 1. However, if fc = 1 or p = fc or p = 0 we obtain 
a degenerate distribution. From here if follows that {qn} is a probability mass function if 
and only if0<p<K<l. □ 


The distribution given by the probability generation fun ction Q with to = 1 is some¬ 
times called modified geometric di stribution ( Phillips! (1978)) or zero-modified geometric 
distribution ( Johnson et al.l ( 2005h L This distribution is obtained as a mixture of a de¬ 
generate distribution and geometric distribution: let 17 be a degenerate random variable 
identically equal to zero, and let T be a geometrically distributed random variable with 
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parameter b G (0,1]. Let q G (0,1) and denote Z = qU + {I — q)V. Then the probability 
generating function of the mixture Z is given as 


Q{z) = q+{l-q) 


bz 


l-{l-b)z' 


We can reparametrize this distribution, by putting 


q = - - — and b 

1 — pn 


and b = 


1 — K 
1 — pK 


with p G (0,1) and k G [0,1). Then the probability generating function takes form (12.1211 
with m = 1. 

The parameter m specifies the lattice of the distribution. We will denote the distribu¬ 
tion with probability generating function Q by G{p, k, m). If m = 1, we will write simply 
Q{p, k). 

Lemma 2.10. The function Q{z) can be decomposed as 



Proof. The decomposition can be verified by computation, as 



□ 


The function Bp(z) is the probability generating function of the Bernoulli distribution. 
In previous Subsection we showed, that Bp generates a commutative semigroup. Using 
the decomposition (12.Kill it is easy to see that the semigroup Q is commutative, as 


Qpi = S ^oBp^oSoS ^oBp^oS{z) 

= S~^ o Bp^ o Bp 2 o S{z) 


and we already showed that Bp.^ o Bp.^{z) = Bp^p.^{z). 

If we choose m = 1 and At = 0 the modified geometric distribution reduces to the 
Bernoulli distribution. We can modify the operator Q and consider two-sided thinning 
operator of geometric type. This can be done by considering 


TZ{z) = S ^oBpoS^^\z) where S^^\z) = qS{z) + {I — q)S{z g G [0,1] 


instead of Q{z). We will denote two-sided modified geometric distribution by 2C/(p, k, q, m). 
We see that Q is obtained from TZ by considering q — 1. 

We will study discrete stable distributions with G thinning operator (of geometric 
type) more into details in Sections 3-5. It will be shown there that this choice of thinning 
operator in Definition 12.11 leads to a distribution with probability generating function 
given by 


(2.14) 



A > 0, 7 G (0,1], K G [0,1), m G N. 
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Thinning operator of Chebyshev type. Let us consider a function of the fol¬ 
lowing form 


(2.15) 


Q(^) 


(1+^) 


where p G (0,1) and b G (—1,1) and Tp{x) = cos (parccosa;). 


Remark 2.11. The function T„ for n G N is called Chebyshev polynomial. It belongs 
to the class of orthogonal poly nomial s . The re is an extensive literature about Chebyshev 
polynomials, see for example iRivlinI (1197411 . Chebyshev polynomials are commutative, 
Tn o Tm(x) = TmO Tn{x); they have the nesting property, Tn o T^ix) = Tmnix). This 
holds true also for Tp{x) with p G (0,1), defined as Tp{x) = cos (parccosx). However, in 
this case Tp is not a polynomial any more. 


The function Q{z) can be decomposed in the following way: 
(2.16) 

Q{z) = R~^ o Tp circR{z), where R{z) = ^ 

z — (i + bjz 


2(fe + y) 
(l-t 5 )(l-t?/)' 


Lemma 2.12. The function Q{z) is a probability generating function. 


Proof. Let us consider only the case of 6 = 0 and p='i,nGN, n>2. Then we can 
rewrite the function Q{z) in the following form 


Q(z) = 


2 cos arccos 5^) 

1 -I- cos I - arccos 7^ ) 

Y n 2 — z I 


Using the exponential and logarithmic forms of cos and arccos functions cos(x) = (e‘^ -I- 
e~^^)l 2 and arccos(x) = f + ilog (ix -I- \/l — x^) we can rewrite cos(i arccos j/) into the 
following form 


cos(i arccos y) 


1 

2 


-log (iy+\/i-y^) 


_l_ g-iilr+i°g(iy+\/i-y^ 



. / _\ —Xjn . / _\ 1/n 

= [iy + \/l - + 56“*^ Y?/ + 

_ i e‘S +(iy+ x/l-y^)^/" 

^ (iy -I- \J\ — y2)i/" 


Hence Q{z) simplifies into (we use substitution y = 5 ^) 

cw = 2 + (i» + %/iT g)Ll_ 

gi^ _|_ (jy _|_ y^l — y2^1/ra 

^ 0 l + (y-i\/l-y^)^/" 
l + iy- ii/l - 
2 

“ 1-H ^ 


So for z G (0,1] we have 
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We have to show that Q{z) is a real function of z. Let x = 2 ^, y = —2and 
u = X + iy = r(cos (f> + i sin (j)). Then using Moivre’s formula 


2-z 


-2i 


. \/l — 


2-z 


Ijn 


’ -2i'^ 


— 1 /n 


^2 — z 2 — z 
= r^^"'(cos(0/n) + isin((/)/n)) + r“^/"(cos((/)/n) — isin((^/n)). 


This number is real if and only if r = 1. But 

. ^ + 4(1 — 2:) 

r = \\x + iy\\ = ^yx^+ y^ = — - -— = 1. 

(2-z)2 

We conclude that for 2 : G (0,1] the function Q(z) is real valued. Moreover Q(l) = 1. 
To complete the proof we need to show that Q{z) is a power series with nonnegative 
coefficients expressing probabilities. 

We denote Q related to the parameter p by Qp(z). The inverse function of Qp{z) is 


2Tr, 


Qp iy) = 


(*) 


1 + t; 


(*)■ 


This follows from the decomposition (12.1611 . Qp{z) = R~^ oTpO R{z), where R{z) = 
and from the fact that the inverse function of Tp{x) is T„(x). This can be verified easily 
from the definition Tp{x) = cos (parccosx). For n G N is T„ the Chebyshev polynomial. 

Consi der first the simple case of n = 2. We know that T 2 {x) = 2x^ — 1 (see, for 
example, RivlinI ( 1974lL l. Therefore 


Qp\y) = 1 +-j- 


y y" 


We may inverse this function again to obtain 


Qpiz) = Qi/2iz) = 


-2 + 2 ^ 2 ^ 
1 - 2 


for 2 < 1. The power series expansion is now easy to obtain 


Qi/2iz) — ^ 


m—0 


V2 

2m,+ l 


(-1)^ 


2 

m + 1 


2F1 (1, i + TO, 2 + TO, \) 2" 


It can be verified that the coefficients 

72 


Pm — 


2rn.+ l 


(-1)^ 


2 

m + 1 


2 F 1 (l,i+TO, 2+TO, i) 


are all positive as is positive for to even and negative for to odd and the hyperge¬ 

ometric function 2 ^ 1 ( 1 , ^ + to, 2 + to, i) is always positive for to > 0. Therefore 21 / 2 ( 2 :) 
is a probability generating function. 


Now we will show by induction that Qpiz) is a probability generating function for all 
p of the form p = 1/2^, with fc G N. We already showed that it is true for p = ^ Let 
us assume Qpiz) is a probability generating function for p = ^, A: > 1. Because of the 
nesting property of Tp we have Tp /2 = Tp oTi/ 2 -, therefore we may write 
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2p/2(-z) = R ^ O Tp /2 O R{z) = R ^ oTpO T 1/2 o R{z) = 

= R~^ oTpO Ro R~^ o Ti/2 o R{z) 

= Qp O Qi/2(-Z) 

By induction assumption Qp{z) is a probability generating function, as well as 21 / 2 ( 2 )- 
The composition of two probability generating function is a probability generating func¬ 
tion itself, therefore we conclude that Qp/ 2 iz) is probability generating function. □ 

We denote the probability distribution given by the probability generating function 
(b^-151) by T(p,6). 

Proposition 2.13. Let eT{p,b). ThenF,e = p^. 


Proof. We compute the expectation of e using the property of probability generating 
functions as Ee = 2^(1)- By deriving 2 ( 2 ) we obtain 

^ 2{1-b)^Tp{uiz)) 

> {l + b){l + Tp{u{zW^ 

where 


u{z) 

^Tp{u{z)) 


u'{z) 


(1 -I- b)z — 2b 
2-{l + b)z' 

^Tp{u)u'{z), 

2{l + b){l-b) 

(2-(l + 6)z)2- 


Using the relation b etween Chebyshev polynomials of the first and second kind (see 


Erdelvi et al.l (ll958alB we obtain 


_d 

du 


Tp{u) = pUp-i{u) = p— 


sin(p arccosit) 


sin(arccos u) 

Putting all together and setting z = 1, u = u(l) = 1 we obtain 

4(l-6)p2Utfc 

n'(i) = J ^ 2 
’ 4(1 + 6) ^ 


□ 


The semigroup Q generated by probability generating functions of this form is com¬ 
mutative. From the decomposition (12.1611 follows tat 

Qpi {Qp 2 i^)) = oTp^o Ro R~^ o Tp^ o R{z) 

= R~^ o Tpj o Tp 2 o R{z). 


But 


Tpi o Tp^ (x) = cos {pi arccos (cos {p 2 arccos a;))) 
= cos {piP 2 arccos x) 

= Tp^ °Tp^{x). 
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We will study discrete stable distributions with Chebyshev type (T) thinning operator 
more into details in Section [S] It will be shown there that this choice of thinning operator 
in Definition [Q leads to a distribution with probability generating function given by 

(2.17) P(z) = exp |-A ^arccosyi^^^p^^ 7 S (0,2], A > 0, 5 e (-1,1). 

2.2 On second definition of discrete stable distributions 

In this Subsection we give a definition of discrete stability that generalizes the second defi¬ 
nition of strict stability (EH) for discrete random variables. The constant An in (EH) takes 
form An = for some 0 < a < 2. Hence the product AnX is generally not integer¬ 
valued and we have to find a different normalization. Compared to the normalization 
used in previous Subsection we need a “portlying” normalization rather than thinning, 
therefore we will look for distributions with expected value bigger than 1 . 

Definition 2.14. Let X, Xi,X 2 , ■ ■ ■, Xn, ■ ■ ■ denote a sequence of independent and iden¬ 
tically distributed non-negative integer-valued random variables. Assume that for every 
n G N there exists a constant Pn > 0 such that 

n X 

(2.18) X{pn) = '^Xi, where Xjpn) = y^£j(pn), 

i=i i=i 

and ej{pn) are i.i.d. non-negative integer-valued random variables. Then we say that X 
is positive discrete stable random variable in the second sense. 

Let us denote the probability generating functions of the random variables X and s(pn) 
by 7^(z) = E[z^] and Qp„{z) = respectively. There is an equivalent definition 

of positive discrete stability in the second sense in terms of those probability generating 
functions. 

Proposition 2.15. A random variable X is positive discrete stable in the second sense 
if and only if for all n gN there exists a constant > 0 such that 


(2.19) 


r{QpAz)) = 'P-iz)- 


Proof. It follows from the definition (12.1811 that X is positive discrete stable in the second 
sense if and only if 

The probability generating function of X can be computed in the same way as in Propo¬ 
sition EH We obtain 


00 

rj^{z) = Y,F{x = k) (e 

k=0 


ySlipA 


Hence X is positive discrete stable in the second sense if and only if its probability 
generating function satisfy the relation 

( 2 . 20 ) v{Qp^{z))=V'^{z). 


□ 

Further denote by Q a semigroup generated by the family of probability generating 
functions {Q{z) = Qp„{z),n G N} with operation of superposition. We show that the 
semigroup Q must be commutative. 
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Theorem 2.16. Let X he positive discrete stable random variable in the second sense. 
Then the semigroup Q must be commutative. 

Proof. Let us denote G{z) = log'P{z). Then (12.191) is equivalent to 


(2.21) nG{z) = G{Qp„{z)), n G N. 

Let G{z) be a solution of (12.211) . Then for all n G N it must hold 

QpAz) = G-^ {nG{z)). 

It follows from here that for all ni,n 2 G N 

(Qp„, (z)) = G-i (niG (G’^ (n^Giz)))) 

= G~^ (nin2G(z)) 

— Qpn2 (2pni (^)) J 

which means that Q is commutative. □ 

In the following Subsections we introduce some examples of commutative semigroups 
Q leading to several possible distributions that are discrete stable in the second sense. 


Degenerate portlying operator. Assume that the probability generating function 
Q(z) = z", i.e. the portlying distribution is a degenerate one taking only one value n. It is 
obvious that the semigroup Q is then commutative. This choice of Q leads to a distribution 
with probability generating function V{z) = z, i.e. a degenerate distribution localized at 
point 1. We are dealing with a simple summation n = 


Geometric portlying operator. Let us consider now geometric distribution with 
parameter p G (0,1) with probability generating function 


Q(z) = 


pz 


1 - (1 -p)z' 

Such distribution generates a commutative semigroup i 

P 1 P 2 Z 


as 


Qpi (Qp2 (z)) — 


P 1 P 2 Z 


1 - (1 - p2)z - (1 - Pl)p2Z 

— Qp2 (Qpi (z)) ■ 


1 - Z + P 1 P 2 Z 


Proposition 2.17. Let X be an integer-valued random variable with probability generat¬ 
ing function 

P(z) = exp|-A(^l-i^ |. 

Then X is positive discrete stable in the second sense. 


Proof Let Q{z) = and set p so that p z = n. Then 


logP(Q(z)) = -A 1- 


1 - {l-p)z 
pz 

1X7 


= -A 


- 1 + (1 - p)i 
pz 


= — Ap '*'1 -) = nlogP(z). 


Hence by Proposition 12.151 the random variable X is positive discrete stable in the second 
sense. □ 


It is important to note that the probability generating function V{z) defines a non¬ 
positive integer-valued random variable. 
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Portlying operator of Chebyshev type. Consider a probability generating func¬ 
tion 


( 2 . 22 ) 


Q(z) = 




n e N, 


where Tn{x) is the Chebyshev polynomial, Tn{x) = cos(narccosa:). iKlebanov et al. ( 2012 1 
showed that the function Q{z) = Qn{z) is indeed a probability generating function of 
a random variable with values in N. The semigroup Q generated by the family {Q„(z), n G 
N} is commutative. We have Q{z) = R~^ o S o R{z), where R{z) = ^ and S{x) = T„(a;). 
Hence 


Qni {Qn2 (z)) = R ^ oTn^O Ro R ^ O O R{z) = R ^ oTn^oTn^O R{z) 

= R-^ oTn^o Tni O R{z) = Qn 2 (Qm (z)), 
because Chebyshev polynomials are commutative. 

Theorem 2.18. Consider the following function 

(2.23) Viz) = ( ^ ~ ~ ^ j : MgN. 

Then V is a prohahility generating function of a random variable on N. Moreover if X 
is an integer-valued random variable with probability generating function V then X is 
positive discrete stable in the second sense. 


Proof. Let us show first that V{z) is a probability generating function. We will consider 
only the case M = 1. For M > 1 the result will follow as V{z) = 'Pf^{z), where 
Vi{z) = and integer power of a probability generating function is a probability 

generating function of a sum of i.i.d. random variables. It is obvious that P(l) = 1. We 
can write Vi^z) as 


r{z) 




The coefficients of the series are all positive, because the binomial coefficient (|) involves 
(/c — 1) negative factors. 

Now let us show that X is positive discrete stable in the second sense. Let Q{z) be as 
in (12.2211 . Then 


P(Q(V) = (i) - (1) - 1 . 

We can use the explicit expression of Chebyshev polynomial to obtain 

_ (1 + + (1 - 

'^\z) 2z^ 


and 




16 

























From here we see that 


V{Q{z)) = V{zT. 

Hence by Proposition 12.151 the random variable X is positive discrete stable in the second 


sense. 

In the proof of the theorem we showed that 


□ 




k=l 


-If 2 
k 


^2k-l 


SO the probabilities P(X = k) are given as (—1)^ ^(|) for all odd k > 0 and 0 otherwise. 

Remark 2.19. The probabil ity distribution with generating function (12.21111 for M = 
1 is known (see ( Felleil 19681 §XI.3)) as a distribution of the first passage time of a 
random walk through +1. Let us consider a sequence of Bernoulli trials Xi, X 2 , ■ ■ ■ 
with probability p = 1/2, i.e. P(Xi = 1) = 1 — P(Xi = —1) = 1/2 and denote Sn = 
Xi + X 2 + ■ ■ ■ + Xn, So = 0. Then the random walk S'„ passes through +1 for the first 
time at time m if 

<0, S^ = l. 

The probability of this event is given by the probability generating function (12.231) . 

In continuous case we have a similar result. The first passage time of a Brownian 
motion through a level a > 0 has Levy distribution, a special case of stable distribution 
with a = 1/2. 

The discrete stable distribution with probability generating function (12.231) with M = 1 
can be considered a discrete analogy of Levy distribution as is shown in the following 
Theorem. 

Theorem 2.20. Discrete stable random distribution in the second sense with probability 
generating function 


V{z) = 


1 - VT^. 


belongs to the domain of normal attraction of stable distribution S(^, 1,1,0), i.e. Levy 
distribution. 

Proof. Let Xi, X 2 ,..., Xn be i.i.d. positive discrete stable random variables in the second 
sense with probability generating function V(z). The characteristic function of X\ is equal 
to f{t)=r{D*). Denote 

1 


Sn=—y,X,. 

i=l 

Then the characteristic function of Sn is equal to 

fn{t) = P {t/n^) —)> exp |-\/2(-it)^/^|, as n-)■ 00 . 
Moreover (—(1 — isgn(t)). 


□ 


Consider now a slightly different setting with portlying operator with probability gen¬ 
erating function 


(2.24) 


Q{z) = 


1 


llm 


Tn 


, n, m S N. 


As was noted in iKlebanov et al.l (l2012f) . Q is a probability generating function of a random 
variable with values in mN. 
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Theorem 2.21. Let X be an integer-valued random variable with probability generating 
function 


V{z) = 


1 — \/l — z 


2m 


TO S N. 


Then X is positive discrete stable in the second sense. 
Proof. We have 





and using results from the proof of Theorem 12.181 


r{Q{z)) = 


1 - \/l - 


□ 


2.3 On third definition of discrete stable distributions 

In this Section we give a definition of discrete stability that generalizes the third definition 
of strict stability (ESI for discrete random variables. As it turns out, this definition is 
a combination of the two previous definitions. 

Definition 2.22. Let X,Xi and X 2 be independent and identically distributed non¬ 
negative integer-valued random variables. Assume that for any positive numbers pi and 
P 2 there exists a positive number p such that 

(2.25) X{p) = Xi{pi) X 2 {p 2 ), where X{p) ='^ej{p) 

3 = 1 

and ej{p) are i.i.d. non-negative integer-valued random variables. Then we say that X is 
positive discrete stable random variable in the third sense. 

Let us denote the probability generating functions of the random variables X and e{p) 
by 'P{z) = E[z'^] and Qp{z) = E[z®^p^] respectively. Let us again denote the semigroup 
generated by {Qp,p G A} with operation of superposition by Q. There is an equivalent 
definition of positive discrete stability in the third sense in terms of those probability 
generating functions, following directly from the Definition. 

Proposition 2.23. A random variable X is positive discrete stable in the third sense if 
and only if for any positive numbers pi and p 2 there exists a positive number p such that 

(2.26) V{Qp{z)) = V{QpA^))V{QpAz))- 

We can show that every random variable positive discrete stable in the first sense is 
also positive discrete stable in the third sense. 

Theorem 2.24. Let X be positive discrete stable in the first sense. Then X is positive 
discrete stable in the third sense. Moreover (12.2511 holds with 

= pI + pI- 


18 











Proof. Let X be positive discrete stable in the first sense, and let Xi,X 2 ,... be inde¬ 
pendent copies of X. Then the semigroup Q is commutative, p S A = (0,1) and for any 
n > 2 there exists a constant Pn G ( 0 , 1 ) such that 

n 

From Theorem 12.61 we know that pn = Let pi,p 2 ^ A- Then for all ni,n 2 > 2 

ni ni+n2 

Pi(z)Xi+p2QX2 = 'Y^PiPmG)Xi+ ^ piPn^QXj. 

i—1 

If pJ^P 2 cire rational, then we can find n 2 ,p such that 


PlPni —PPni+n2; 
P2Pn2 ~ PPni~\-n2-} 


or equivalently 


■y y 

Pi =p^ - 


pI =p^- 


ni + n2' 
n2 

ni + 712 


But then, with n = ni + n 2 

n 

Pi O + P2 O ^2 = ^ PPn Q Xi = pQ X. 


Moreover pi,p 2 ,p satisfy the relationship p'l+p'l = p^. By continuity argument it follows 
that (12.2511 hold for any choice of pi,P 2 with p such that pf + pif = p'^. □ 


Under some additional conditions we may show that the opposite statement holds true 
as well. 

Theorem 2.25. Let X be positive discrete stable in the third sense and assume that the 
semigroup Q is commutative, A = (0,1) and that there exists a constant 7 > 0 such that 

P^ =p'l +P2- 

Then X is positive discrete stable in the first sense. 


Proof. We may show this by induction. Because X is positive discrete stable in the third 
sense, we have for pi = P 2 = 2~^l^ that 

X = Xr{p2)PX2{p2). 

Let n > 2 and let us assume that 

n 

X = '^Xi{pn), with = 

Denote Y = ^i{Pn) and let p = (j . Because X is positive discrete stable in 
the third sense, Y = X and p'’' + pf^_^_l = 1, we have 

A = y(p) + A„+i(p„+i). 
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The probability generating function of the right-hand side is 


Vy{Qp{z))V{Qp^^,{z)) = 

= P"(Qp„pW)7’(Qp„^,(^)) 

= 7^"+'(Qp„+.W), 

because Q is commutative and pp„ = Pn+i- Therefore 

n+1 

i=l 

□ 


Binomial thinning operator. Let us consider the case of the binomial thinning 
operator with probability generating function Q(z) = {l—p)+pz. Then a random variable 
X with probability generating function V{z) = exp {—A(1 — 2 )^} is positive discrete stable 
in the third sense, as (12.261) holds if 

=p7 +P2- 


Modified geometric thinning operator. We can verify that the positive discrete 
stable random variable in the first sense with modified geometric thinning operator is also 
positive discrete stable in the third sense. Let X be a positive discrete random variable 
in the first sense with probability generating function V{z) = exp ^ Thsn 

P(C,W) = exp|-Ap-'(^^) }■ 


Thus (12.261) holds if 


= pI + pI- 


Chebyshev thinning operator. In the same manner we see that a positive discrete 
stable random variable in the first sense with Chebyshev thinning operator X is positive 
discrete stable in the third sense. Let V be as in (12.171) and Q as in (12.151) . We have 

P(Qp(z)) = [r{z)f. 

Therefore again (|2.26l) holds if 

p'^ =p'l +pI- 


Chebyshev portlying operator. Now let’s look at an example with Chebyshev 
portlying operator with probability generating function Qn{z) = l/Tn{l/z). Then a ran¬ 
dom variable X with probability generating function V{z) = (l — \/l — 2 ^) /z is positive 
discrete stable in the third sense, as (12.261) holds if 


n = ni -I- 712. 


3 Properties of positive discrete stable random vari¬ 
ables 

Distributions, that are discrete stable in the first sense, form the widest and most in¬ 
teresting class of distributions, and in the following Sections we study properties of the 
distributions with thinning operator of geometric type. 
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To remind the definition, a non-negative integer-valued random variable X is said to 
be positive discrete stable in the first sense, if 

n Xj 

(3.1) where 

i=i i=l 

where Xi,X 2 ,... are independent copies of X and el are i.i.d. non-negative integer¬ 
valued random variables. Throughout this Section we will assume that the random vari¬ 
ables come from modified geometric distribution Q{p, K,m) with probability generat¬ 
ing function Q of the form 


(3.2) Q(z) = 


{I — p) + {p— k)z^ 
(1 — pn) — k(1 — p)z’'‘ 


{0 < K < 1, 0 < p < 1, 771 = 1} 

or 

{0<p<K<l, ?n G N, TO > 1}. 


We remind that Q(z) can be decomposed as Q{z) = S ^ o Bp o S{z), where Bp{z) = 


pz + (1 — p) and 


*S'(^) = 


(1 - k)z" 
1 — 


S-\y) = 


{1 — k) + Ky 


Theorem 3.1. A non-negative integer-valued random variable X is positive discrete stable 
with Q thinning operator if and only if Q takes form (13.21) and the probability generating 
function V{z) = Ez^ is given as 


(3.3) 'P{z) = exp < —A 


1 - z^ 

1 — KZ" 


m \ 7 


with 7 G (0,1], A > 0, K G [0,1), 771 G N. 


Proof. Let h(z) = logP(z). From Proposition 12.21 it follows that X is positive discrete 
stable if and only if h{z) = nh{Q{z)) for all n. Set 

h{z) = -\(^^^ 


KZ'' 


and select 7 such that l/p'>' = n. We see that 

1 - z™ , (1 - k)z^‘ 


= 1 - 


1 — KZ'" 1 — TtZ’^ 

Therefore, using the decomposition of Q(z), 


= 1-S{z). 


nh{Q{z)) = —Xn (1 — S'(Q(z)))’^ = —Att (1 — i?p(S'(z)))" 
= —Xn (p — pS{z)y = —Xnp'^{l — S{z)y 

1 _ z"? \ ^ 


= -A 


1 — KZ’^ 


= h{z). 


□ 


The parameter m determines the size of the lattice of the distribution. We will denote 
positive discrete stable random variable (and associated distribution) by PDS’"( 7 , A, k). 
In the case when m is omitted we will understand that ? 7 i = 1. If moreover n is omitted, 
we will understand that k = 0 , in which case the disc rete stable distribution redu ces to 
the discrete stable distribution as it was introduced in Steutel and van Harn ( 19791) . 
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The characteristic function is given as 


f{t) = exp 


f 1- 
\ 1 - 



The case of 7 = 1 is a special one as it leads to a distribution with finite variance and 
exponential tails. As a simple corollary we obtain Poisson distribution by taking k = 0 
and 7 = 1 . 


3.1 Characterizations 

In this Subsection we present several characterizations of positive discrete stable random 
variables. 

Theorem 3.2. Let 7 G (0,1) be a given parameter. Let X, Xi, X 2 , ■ ■ ■ be i.i.d. non- 
negative integer-valued random variables and Y be a non-negative integer-valued random 
variable, independent of the sequence Xi, X 2 , ■ ■. ■ Then X is positive discrete stable 
PDS( 7 , A) random variable if and only if 

r X 

(3.4) A = ^ 0 Aj, where pQ X = £i{p) 

j=i i=i 

and Si{p) are i.i.d. Bernoulli random variables with probability generating function Qp{z) = 
1 — p + pz. 


Proof. First let us show that if A is PDS( 7 , A) then it has the representation (13.41) . Let 
V{z) be the probability generating function of A. The probability generating function of 
the right-hand side of (EH) can be computed in the following way. 


E 





= E 


E 



y-^^->qx,^Y 


^HiQY-U.iz))] 


= E[exp{-Ay (1 - QY-Y.{z)f}\ = E [exp{-Aryi (1 - z)^}] 
= exp{-\{l-zf} = V{z). 


The proof of the inverse statement is more complicated and relies on the method if 
intensively monotone operators. The condition (El can be translated into the form of 
probability generating functions as 

00 k 

(3.5) V{z) = Y,ny = k)X{V{Q^-uAz))- 

k—0 j—1 

Put G{z) = \ogV{z) and h{z) = G(z)/(1 — z)'*'. Then we can rewrite (13.51) as 

00 k 

(3.6) h{z) = (1 - z)-^ Y. = ^) E (1 - Qk-u-, {z)f h (Qfc-v. {z)) 

k—0 j — 1 

00 

= (1 - z)-"^ ^ P(y = k){l- zf h (Qfc-v. (z)) 

k^O 

00 

= YnY = k)h{Q^-u.{z)). 

k—0 

Let A be an operator acting on 5 G G[0,1] such that 

(Ag)(z) = I " k)g(Q,-.,.(z )), z < 
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We can verify that A is an intensively monotone operator (see iKakosvan et al. ( 1984 )) 
and t hat Aff g CjO,!] . It is clear that Aa = a for all constant functions a. It follows 
from ( Kakosvan et al. . 19841 Theorem 1.1.2) that the only solution of (IT (ill is identically 
equal to a constant. Hence h{z) = —X and 


Viz) = exp{—A (1 — z^} . 


□ 


Theorem 3.3. Let 7 , 7 ^ G (Oj 1] assume that 7 ' < 7 . Let S.y be a 'y-stable random 
variable with Laplace transform exp{—u'*'}. Then 

PDS( 7 ', a, k) = PDS ( 777 , k) . 

Proof. The characteristic function of the right-hand side can be computed as 

\ y/'l' 


E 


exp |itPDS ^777) A^^'^S.y, n'j | 


)}] 

= E 

exp \ 



L 


= E[exp {itPDS (Y, A, k)} ]. 



□ 


The following Corollary can be applied for simulations of positive discrete stable ran¬ 
dom variables. 

Corollary 3.4. Let Y,Yi,Y 2 ,... be a sequence of i.i.d. random variables with geometric 
distribution, P(y = n) = (1 — n > 1. Let N be a random variable, independent of 

the sequence Yi, Y 2 ,.. ■, with Poisson distribution with random intensity where 

S.y is a "f-stable random variable with Laplace transform exp{— m"''}. Then 

N 

1=1 

has the same distribution as a positive discrete stable random variable PDS( 7 , A, k). 

Proof. Let X = Then X is a compound Poisson random variable with random 

intensity A^/'^'S.y and jumps Yi,Y 2 ,... with characteristic function 


9{t) = 


(1 — K)e 

1 — KC*' 


The characteristic function of a compound Poisson random variable with intensity r 
and characteristic function of jumps h{t) is exp{— t(1 — h{t))}. Therefore X is in fact 
PDS(l,Ai/^S.^, k). We thus obtain the result from the previous Theorem 13.31 with 7 ' = 
7 . □ 


3.2 Moments 

Theorem 3.5. Let X be PDS( 7 , A, k) random variable with 7 = 1 and k > 0. Then the 
n-th factorial moment can be computed using the following formula 


(3.7) 


E [(X)„] = 


(l-«)^ 


n —1 

s=0 


1 


7+1)! 


- lA A*+1 


^s+l 


23 


























Proof. Let V{z) be the probability generating function of X. The n-th factorial moment 
of discrete random variable can be computed as the value of the n-th derivative of the 
probability generating function at point 1, i.e. 

e™= _ . 


Since P{z) = exp{(/( 2 :)}, with 


g(z) = -A 1 - (1 - k) 


1 — Kzy 

we compute the n-th derivative using the Bruno’s formula ( Faa di Brunol ( 1857t) ~) 


d" 

dz^ ^ ’ 


Z = 1 






where Bn k{xi ^..., Xn-k^i) is the Bell’s polynomial, 
(3.8) 


Bn,k ■ ■ ■ : fc+l ) — ^ ^ 


n\ 


Zl,...,Zn-fc + l 


il\i2\ . . . in-k+l 


:(!()•■ (t)" 


^n—/c+1 

{n — k + 1)! 


-fc + 1 


where we sum over all possible combinations such that ii + 2 i 2 ■ ■ ■ + {n — k+ l)in-k+i = n 
and ii + 12 - ■ ■ + in-k+i = k. By differentiating the function g{z) we obtain 


#(l) = z!A- 


..i-l 


(1 — k)® 

Plugging that into the Bell’s polynomial we obtain 


T(^)' 




zi!z 2! .. .z„-fc+i! V j' J 


= E 


n\ 


zi!z 2! .. .z„-fe+i! \(1 -k) 


J- 

n—k-\-l 

n 


A/c-l 


n! A^n" 

K^(l — k)’ 
k 


Zi!z 2! • . .Zn-fe+l! Kf{l - k)” 

-B„,fe(l!,2!,...,(n-fc+l)!) 


Xk 


X{1 — k)" \k j \k — 1 


n\ fn — 1 


(n — A:)!. 


Hence the n-th factorial moment is 

n 

e[{xu = Y^ 


k 


Xk 


k=l 


n\ fn — 1 


X(1 — k)" \k j \k — 1 


^ A^ n! /n — 1 


(n — fc)! 


(1 — kX ^ X k\\k — l) 

^ ’ k=l ^ ^ 

The result follows from here by setting s = fc — 1. 


□ 
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3.3 Probabilities 


In the next Theorem we show connection between the probabilities of a positive discrete 
stable random variable and moments of a tempered stable random variable. 

Theorem 3.6. Let X be a PDS( 7 , A) random variable with 7 < 1. Let Y be a tempered 
stable random variable with characteristic function /y(t) = exp{—(A^/'*' — it)''' + A}. Then 
we can write the probabilities P(X = k) as 

\^h 

F(X = k) = 

k\ 

Before we proceed to the proof of the Theorem, we state a simple Lemma. 

Lemma 3.7. Let be j-stable random variable with Laplace transform L{u) = = 

exp{—rt'*'} and density function p{x). Let 6 > Q. Let Y be a random variable with density 
function 

prix) = e~^^p{x)/L{9). 

Then Y is a tempered stable random variable with characteristic function 

f{t)=eM-{0-ity + 9^}. 

Proof. We may compute the characteristic function of Y as follows: 

pOO pCO 

/y(t) = Ee‘‘^= / e**"py(x)dx= / e“"e-%(:r)/L(0)da: 

Jo Jo 

poo 

= e® / exp {—(0 — it)a:}p(x)da: 

Jo 

= e^~'L{9 — it) = exp{—(0 — it)'’' + 9^} . 


□ 


Now we can prove the Theorem. 


Proof of Theorem \y.b\ It follows from Theorem l3.3l that a positive discrete stable random 
variable PDS( 7 , A) is a Poisson random variable with random intensity Sj, where 
is a 7 -stable random variable with Laplace transform L{u) = exp{—u'’'} and density 
function p{x). Therefore the probabilities P(X = k) can be computed as 


F{X = k) 


,(Ai/'’'s)* 


fc! 


■p(s)ds 


\k/l 

~kr 


L(A^/^) 




’p(s)/L(A^/'’')ds. 


But e ^ ^^®p(s)/L(A^/'’') is a density function of a tempered stable random variable Y 
with characteristic function /(t) = exp{ —(A^/'’' — it)'*' -|- A}. Therefore 


F{X = k) 


\k/-f pco 

^L(Ai/^) sV(s)ds 
xfc /7 

L(A^/T')Ey'=. 


□ 
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Theorem 3.8. Let X he a PDS( 7 , A, k) random variable with 7=1 and k > 0. Then 
the probability P(AC = m) for m > 1 can be computed using the following formula 


m—1 


(3.9) 


'{X = m) = e 


/m-1 


s =0 


(s + 1)! 




s+l 


Proof. We compute the probabilities by expanding the probability generating function 
into power series. 


Viz) = exp —A ( 1 — (1 — k) 


1 — KZ 




A^ 




n—1 
00 00 


(1 - kzY 


X X 




n—1j—0 
00 00 


= e 


-A I ^-A 




n—1 m—n 
00 m 


= e 


-A I ^-A 


EE^(i-r 


m—1 n—1 
00 m— 1 


n + j — 1 
3 

m—1 
m — n 

m—1 
n — 1 


yn +3 


= e 


-A I „-A 






^ (s+ 1)! 

m=l s =0 ' ' 


The probabilities P(X = m) are obtained from this results as the coefficients of the 
probability generating function by z™, as P(z) = ~ m)z"*. □ 

Corollary 3.9. Let X be PDS(7, X, k) random variable with 7 = 1 and k > 0. Then the 
probability P(X = m) for m> 1 can be expressed in the following ways 


(X = m) = e-^A(l - k)k^-^ iFi[l-m,2 


1 , 


-A 


and 


P(X = m) = e-^A(l - k)k^-^ 

m 


P-l. 


A , 


where iFi(a, 6,z) is the Kummer confluent hypergeometric function and L^fi\z) is the 
generalized Laguerre polynomial. 

Proof. The first assertion follows directly from (IXnil . The second assertion follows from 
the relation betwe en Laguerre polynom ial and Kummer confluent hypergeometric function 
(see for example ( Erdelvi et ah . 1953bL pp. 268)), stating that 


LifHz) = 


n a 
n 


iFi{-n,a + l,z). 


□ 


3.4 Continuous analogies 

Let us consider a random variable X“ = aX, with X ^ PDS( 7 , A, k) and a > 0. Then 
X“ takes values in oNq = {0, a, 2a, - ■ ■}. We study the limit behaviour of X“ as a —>■ 0 
with K —>■ 1 . 
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Theorem 3.10. Let X be a positive discrete stable random variable with parameters 7, 
A and n and let = aX with a > 0. Let k = 1 — ac. Then 


r(() = exp|-A(i-^) } ^ ,,(()= exp I-A }, 

Proof. The limit characteristic function can be computed in a straightforward way. We 
have 

1 — e‘“* 1 — —iat 

- -r-= --r--as a —>■ 0. 

1 — 1 — e‘“* + ace‘“* —lat + ace'“* 

Hence we have 

.p(0=limexp{-A(—}=exp{-A(^^) }. 


□ 


Next we show that discrete stable distribution on Nq can be considered a discrete 
analogy of stable distribution with index of stability 0 = 7 and skewness parameter 
[5 = 1 . 

Theorem 3.11. Let X be a positive discrete stable random variable with parameters 7, 
A and k and let = aX with a > 0. Let A = 5 /a'’'. Then 


nt) = exp 



/ _ giat 

\ 1 — Ke '“*) j 

— ^{t) = exp|-cr|tp (1 - isign(t) tan (^))} . 


as a ^ 0, 


where a = cos (^) . 


Proof. We have 


Hence 


1 - 


1 - 


1 — Ke'“* 1 — k + k( 1 — e*“*) 


—iat 

(1 — k) — Kiat 


as a —>■ 0. 



5 / -iat V 

-7:;-r- as a —>■ 0 

a'>' \ (1 — k) — mat J 

--— , (—it)'^ as a —>■ 0. 


Finally we notice that 

(—it)T' = |tp(—isign(t))'’' = \tp cos (7r7/2) (1 — isign(t) tan (7r7/2)). 


□ 


3.5 Asymptotic behaviour 

In this Subsection we show that the tails of discrete stable PDS(7, A, k) distribution are 
heavy with tail index 7. 
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Proposition 3.12. The discrete stable distribution PDS( 7 , A, k) belongs to the domain 
of normal attraction of a-stable distribution with characteristic function 


g{t) = exp - 


(1 - k)' 


■ cos (7r7/2) Itp — isign(t) tan 


(?) 


Proof. Let ACi, X 2 ,..., Xn be i.i.d. PDS( 7 , A, k) random variables with characteristic func¬ 
tion 


f{t) = exp|-A 

Let us denote Sn the normalized sum 

X1 + X2 


1 - 

1 — Ate'’ 


it \~> 


Sn = 


X„ 


ph 


Then the characteristic function of Sn is given as 


E [e""-] = /" 


1/7 \ ' 


= exp < —A 


1 - e'*/” 
1 — 


We use the Taylor expansion of exp to obtain 

—it 


logE = —Xn 


(1 — K)n^l^ 




A 


( 1 -At )7 


(—it)'’'(l-I- 0 (n as n —>• 00. 


Hence 


g(f) = hm E = exp ( 

n-Aoo (1 — AtjT^ J 


We can rewrite the exponent using 

(—it)'*' = |t|'>'(—isign(t))'*' = It? cos (7r7/2) (1 — isign(t) tan (7r7/2)). 


□ 


4 Properties of discrete stable random variables 

In this Section we will study more into detail the discrete stable distribution in the limit 
sense with two-sided modified geometric thinning operator, as defined in Section [2.11 To 
remind the definition, an integer-valued random variable X is said to be discrete stable 
in the limit sense, if 

n X- 

(4.1) X = Yrni '^Xifpn), where ^ 

i=l j=l j=l 

where Xi, X 2 , ■ ■ ■ are independent copies of X and s^j \ are i.i.d. integer-valued random 

variables. Throughout this Section we will assume that the random variables 4*\ 4*^ come 
from two-sided modified geometric distribution 20{p, K,m,q) with probability generating 
function TZ. We remind that the probability generating function TZ is given as 

(4.2) TZ{z) = S-^oBpoS^^'>{z), 
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where 


S{z) = 


S-^{y) = 


(1 - k)z^- 

1 — KZ™ 


1 - k(1 - y) 


Bp{z) = l-p + pz, 


and finally 

5(2)(z)=95(z) + (1-9)5(z-1). 

Theorem 4.1. An integer-valued random variable X is discrete stable in the limit sense 
with two-sided modified geometric thinning operator, if and only ifTZ{z) takes form (14.211 
and the probability generating funetion 'P{z) = Ez'^ = J2T=-oo = k)z^ takes form 


I + (1 — kIz"* (1 — /c)z 

(4.3) r{z) = exp -A 1 - - (1 - qy 


1 - KZ" 


1 — KZ 


-A 


1 - 


2 J \ 1 — KZ' 

with 7 s (0,1], A > 0, K G [0,1),/3 G [—1,1], g G [0,1]. 


^ ^^(1-k)z™ ^(1 - k)z-'"^ 


1 — KZ" 


Proof. We have shown in Proposition 12.81 that a random variable X is discrete stable in 
the limit sense if and only if 

V{z)= lim [iPo + lPl(7^(z)) + iP 2 ( 7 ^(l/z))], 

n^oo 

where Pi is the generating function of the sequence {pi,P 2 , • ■ ■ } with pk = P(X = k) and 
V 2 is the generating function of the sequence {gi, ( 72 , ■ ■ ■} with q^ = P(X = —k). Let us 
assume that Vi and P 2 take the following form 


(4.4) 


Vfz)=Vfl)-\, 


with 7 G (0,1]. We notice that 


1 - z"’ 

1 — KZ'"' 


1 - z"* 

1 - KZ" 


1 - z”* 

1 - KZ" 


, ^=l, 2 , 


= 1 - 5 (z). 


This simplifies the computation, as 1 — S fR{z)) = 1— (l — p + pS'(^)(z)) = p (l — ^(^^(z)) 
and similarly for 1 — S{P{l/z)). 

We can now compute the limit 


V{z)= lim [Po + Pl(7^(z)) + P 2 ( 7 ^(l/z))]”. 


Let p = n Then 


V{z) = lim [1 - Ai (1 - S{TZ{z))y - A 2 (1 - S{TZ{l/z))yf 

n—>-oo 


= lim 

n—>-oo . 


1 - AipT- (1 - s^^Hz)y - a 2 p^ (1 - s^^Hi/z)y 

, , (1-k)z- ,(1-k)z-™^^ 

= exp <j -Ai 1 1 - q- -—- (1 - g)- 


—A 2 1 — g 


1 — KZ ™ 
( 1 - k)z™^^ 


1 — KZ 

(izzhr: — 

1 — KZ“™ ^ 1 — KZ' 

By setting A = Ai + A 2 and f = we obtain the desired result. 


□ 
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We will denote discrete stable distribution (and random variable) by DS'"( 7 , /3, A, q, n). 
The parameter m specifies the size of the lattice of the distribution. If we omit m then it 
is understood that m = 1. If k is omitted we will understand that /c = 0. If moreover q is 
omitted we will understand that <7 = 1. In this case the probability generating function 
(lOl) reduces to 


exp —A ( —1(1“ ~ ^ 


1- 


a-i/zr 


which corresponds to the discrete stable distribution introduced in iKlebanov and Slamova 
( 201,lh . In the case of /3 = I and q = 1, the DS( 7 ,1, A, 1, k) random variable correspond 
to positive discrete stable random variable PDS( 7 , A, k). 


Remark 4.2. A discrete stable random variable X ^ DS( 7 , P, A, q, k) is infinitely divisi¬ 
ble, as for all n G N, 


X = Yi+Y 2 + - \-Yn, where ~ DS( 7 ,/3, A/n, q, k), i = 1,..., n. 


For the sake of simplicity we will denote 


(4.5) 

(4.6) 


9{z)= (^1-q- 
h{z) = g{z~^) 


(I-k)z'" ^^{1-k)z-"^''^ 


1 — KZ''' 


= l-(l-q) 


1 — KZ 


(1 — n)z 

1 - KZ" 


- 9- 


(1 — k)z 
1 — KZ~' 


7 


Then the probability generating function of a DS( 7 , /?, A, q, k) random variable can be 
written simply as 


V{z) = exp |-A 9 {z) - A ^(^)| ■ 


4.1 Properties 


Discrete stable distribution shares many interesting properties with stable distributions. 
In this S ubsection we show that analogies of Properties of stable distributions (see, for 
example, Samorodnitsky and Taaaul (1994)) hold also for discrete stable distributions. 


Property 4.3. Let Xi and X 2 be independent random variables with Xi ^ DS( 7 , Pi,Xi, q, k), 
i = 1,2. Then Xi + X 2 ~ DS( 7 , P, A, q, k), with 


A — Ai -|- A2, 


PiXi + P2X2 
Ai -l- A2 


Proof. Using the notation (I4.5II - (I4.6L the probability generating function of Xi, i = 1,2, 
is 

Vi{z) = exp |-A* 9 iz) - X, • 

The probability generating function of Ai -|- A 2 is a product of the single probability 
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generating functions. Therefore 


\ogVx^+xAz) = - Ai ( g{z) - Ai P ) h{z) 


— Aq 


2 

1 + /32 


9{z) - A 2 ( ) h{z) 


(\ , \ f^ + A2^2'i ! ^ 

-A|t±^'lgW-Afl^'| h(z). 


where A = Ai + A 2 and /? = (/3iAi + + A 2 ). □ 

Property 4.4. Let X ~ PDS( 7 , A, k). Let a G (0,1). Then X{a) ^ PDS( 7 , a'^X, k). 
Proof. The probability generating function of X(a) is equal to 

exp{-A(l - S'(Qa(z)))’^} = exp{-Aa'^ (1 - 5'(z))'^} . 


□ 


Property 4.5. Let X ~ DS( 7 , P, A, q, k). Then —X ~ DS( 7 , —/3, A, q, k). 

Proof. This follows from the fact that g(z~^) = h{z), where we use the notation (14.511 - 
(ITO . Then the probability generating function of —X is given as 

= exp |-A ^ ’ 


and this is the probability generating function of DS( 7 , —/3, A, g, k). □ 

Property 4.6. Let X ^ 08 ( 7 , P, A, q, k). Then X is symmetric if and only if g = 1/2 or 

p = o. 

Proof. A discrete random variable is symmetric if and only if 'P{z) = V{z~^). Using the 
notation (H31)-(ITO. and the fact that g{z = h{z), it follows that a discrete stable 
random variable is symmetric if and only if 


-A 



g{z) - A 




Hz) - A Hz)- 


But this holds true if and only if /3 = 0 or g{z) = Hz)- The latter condition is satisfied 
only if g = 1 / 2 . □ 

Property 4.7. Let X be DS( 7 ,/3, A, g, k). Then there exist two i.i.d. random variables 
Yi and Y 2 with common distribution DS( 7 ,1, A, 1, k) such that 


X = Yi 
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Proof. Let Yi, I 2 ~ DS( 7 ,1, A, 1, n). Their probability generating function is 


PW = exp|-A(^) }■ 

Moreover, the probability generating function of Y(p) is obtained in closed form, as Y 
are in fact positive discrete stable random variables. So we have 

Similarly as in the Proof of Theorem 14.11 we can compute that 
ViTlpiz)) = exp{-A(l - S{np{z))y} = expj-Ap'^ (l- —— - (1 - 

I V i — KZ L — KZ 

The probability generating function of the difference Yi(pi) — ^ 2 (^ 2 ) is computed as 

iP(7^p,(z))iP(7^p,(l/z)). 

Putting all together we obtain the desired result. □ 


4.2 Continuous analogies 

Let us consider a random variable X°' = aX, with X ~ DS( 7 , /3, A, q, k) and a > 0. Then 
takes values in aZ = {0, ±a, ±2a, • • • }. We show that the limit distribution of 
is a-stable distribution with index of stability 7 and skewness (3. We study the limit 
behaviour of as a —>■ 0 and g —>■ 1 / 2 . 

Theorem 4.8. Let X be a discrete stable random variable with parameters 7 , /3, A, q 
and K = 0. Let X°‘ = aX with a > 0 and let 2q — 1 k, a as a ^ 0. Then 

nt) = exp |-A (1 - - (1 - - 

■ ^ (^) 

—>■ ip{t) = exp A cos ~ i;5sign(t) tan |, as a —>■ 0 . 

Proof. We may rewrite the characteristic exponent of /“(t) as 

((2g — 1)(—iat))'*'— A (( 2(7 — l)(iat))'’', as a —>■ 0 


iogr(t)«-A 


1 + /3 


and because g (1 + a )/2 we have 



{-itr - A 



{ity. 


To complete the proof it is enough to notice that (—it)'’' = |t|^ (cos— isin^) and 
(it)'’'= |t|'’'(cos+ isin ^). □ 

Remark 4.9. It can be shown that the case of /c > 0 leads to a similar result, the limit 
distribution is again (a-stable with index of stability 7 and skewness /3. 

Proof. □ 
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5 Properties of symmetric discrete stable random vari¬ 
ables 


In the previous Section we studied the general case of discrete stable distribution in the 
limit sense. The symmetric version of such distribution is special case with interesting 
properties and we will therefore study it more into details in this Section. The symmetric 
discrete stable distribution in the limit sense is obtained by considering the symmetric 
two-sided modified geometric thinning operator 2Q{a, k, m). 

Theorem 5.1. A symmetric integer-valued random variable X is symmetric discrete 
stable with symmetric two-sided Q thinning operator if and only if the thinning oper¬ 
ator takes form (14.21) with q = 1/2 and the probability generating function V{z) = 

takes form 


with parameters 7 G (0,1], A > 0, k S [0,1) and m gN. 

Proof. The proof follows from the proof of Theorem 14.11 In the symmetric case we have 
Viiz) = 7 ^ 2 ( 2 ), therefore Ai = A 2 and moreover q = 1/2. The probability generating 
function (lOD thus reduces to (ED- □ 


We will denote symmetric discrete stable distribution (and also random variable) by 
SDS'"( 7 , A, k). In case when m is omitted we will understand that m = 1. If k is 
omitted we will understand that k = 0 , in which case the symmetric discrete stable 
distribution reduces to t he sy mmetric discrete stable distribution as it was introduced in 


Klebanov and Slamoval (|2013f) . 


The characteristic function is given as 


f{t) = exp 





cos (tm) — K 
— 2KCOs(tm) -I- 1 



The case of 7 = 1 is a special one as it leads to a distribution with finite variance and 
exponential tails. 


5.1 Characterizations 

Theorem 5.2. Let 7,7^ G (Oj 1 ] assume that 7' < 7. Let S.y be a "/-stable random 
variable with Laplace transform exp{—u'''}. Then 

SDS(7', a, k ) = SDS (777, A^/t'S^, At) . 

Proof. The proof of the Theorem is done in the same way as the proof of Theorem |X31 CH 

Corollary 5.3. Let Y,Yi,Y 2 ,... be a sequence of i.i.d. random variables with two-sided 
geometric distribution, P(K = ±n) = A(1 — K)td^~^,n > 1. Let N be a random variable, 
independent of the sequence Yi,Y 2 ,..., with Poisson distribution with random intensity 
where S..y is a "/-stable random variable with Laplace transform exp{—rt'*'}. A 
random variable X is symmetric discrete stable SDS(7, A, k ) if and only if 

N 

^'Eb- 
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Proof. Let X = X?- Then X is a compound Poisson random variable with random 

intensity and jumps Li,l 2 , • ■ • with characteristic function 

^ 1 (1 - ^)e‘* , 1 (1 - 

2 1 — Ke‘* 2 1 — Ke“'* 

The characteristic function of a compound Poisson random variable with intensity r 
and characteristic function of jumps h{t) is exp{— t(1 — h{t))}. Therefore X is in fact 
SDS( 1 ,A 1 /t'St, , k) . We thus obtain the result from the previous Theorem 15.21 with 7 ' = 
7 . □ 


5.2 Probabilities 

Theorem 5.4. Let X be SDS( 7 , A) random variable. Then 


00 00 


7 jA A-l 1 / i 


.■nii-ii V ■ / V 2 ' 

In case 7=1 this simplifies to 

F{X = k) =e-^Ik{X), kGZ. 


fc e z. 


where Ik is the modified Bessel funetion of the first kind. 

Proof. The generating function of a discrete random variable taking values in Z is a power 
series, with coefficients equal to probabilities, i.e. 


00 

'Px{z)= Y. nx = k)z'^. 

k— — (yc 


(Note that this series converges only for e < \z\ < 1). Thus expanding (15.11) with k = 0 
into a power series we obtain the probabilities. We use Taylor expansion of exponential 
function, binomial expansion and interchange of sums. 


exp 





EEEM) 

j—Oi—O l—Q 


i 


X 1 


j! 2 


2l-i 




change of notation k = 21 — i and interchange of sums 


00 00 00 

k— — oo i—\k\ j—0 



f ^ \ x 1 k 

j! 2 *^ ■ 


From this the first result follows. Taking 7=1 the first binomial coefficient (fl) turns 0 
for j < i and we have, for k > 0, 


p{x = k) = YYi-^y^' 

i—k j—i 


j\ f i \ X I 
i 


\ 

\ 2 


= e-^^(A/2) 


k-\-2l 


j\ 2 * 

1 


1=0 


r(/ + i)r(^ + fc + 1 ) 


= e-^/fc(A). 


□ 
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5.3 Continuous analogies 

Let us consider a case of random variable X°' = aX, with X ~ SDS( 7 , A, k) and a > 0. 
Then takes values in aZ = {0, ±a, ±2o, • • • }. We study the limit behaviour of as 
a —>■ 0 with K —>■ 1 . 

Theorem 5.5. Let X be a symmetric discrete stable random variable with parameters 7 , 
A and k and let = aX with a > 0. Let k = 1 — ac. Then 


/“(t) = exp<^ -A 1 - (1 - k) 


cos (at) — K 
' — 2ncos{at) + 1 


(f>{t) = exp s — A 




as a —>■ 0 . 


Proof. The limit characteristic function can be computed in a straightforward way. We 
have 


1 - (1 - k) 


cos(at) — K 
— 2KCOs{at) + 1 


= 1 + ac 


1 — cos(at) — 


1 + 


2(1 — ac)(l — cos(at)) + 
acf^ 12 — 


t^ — a cP + 


as a —>■ 0 


Hence we have 


t^ 


I \ j ^ I 1 adf 12 — . , 

*’(*) =‘“.“Pi-^1 +(2_ace+ ed h“Pl-\e + c2 


□ 


Next we show that symmetric discrete stable is a discrete analogy of symmetric stable 
distribution with index of stability a = 27 . 

Theorem 5.6. Let X be a symmetric discrete stable random variable with parameters 7 , 
A and n and let = aX with a > 0. Let A = bjcsf^. Then 


/“(t) = exp<^ -A 1 - (1 - k) 


cos (at) — K 
' — 2 KCOs(at) + 1 


^ I —> ip(t) = exp {—o’ltp'’'} , as a —>■ 0 , 


where a = ^ 
Proof. We have 




— 2 KCOs(at) + 1 


(1 + k) 


— 2 KCOs(at) + 1 
aH^ 


Hence 


-A 1- (1-k)- 


cos(at) — K 
■f — 2 «;cos(at) + 1 


2 (1 — k) 2 + KO^t^ 

b f (1 + k) a^t^ 


as a —>■ 0 


a^T" \ 2 (1 — k)2 + 

a® a ^ 0 . 

2'1 (1 - k ) 27 ' ' 


as a —>■ 0 


□ 
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5.4 Moments 


In this Subsection we give a formula for factorial moments of SDS(1,A, k) distribution 
and show that fractional moments of SDS( 7 , A, k) of non-integer order up to 27 exists. 

Theorem 5.7. Let X be SDS( 7 , A, k) random variable with 7=1 and k> 0. Then the 

n-th factorial moment can be computed using the following formula 

(5.2) 

1 \k 

E = (13^ E (0> 2!(« - 1), 3!(«2 + 1),..., (n - fc + !)!(«"-'= - (-1)"-'=+!)) , 

where Bn,k is the Bell’s polynomial (ES- 

Proof. The proof is analogous to the proof of Theorem 13.51 and therefore is omitted. □ 

Theorem 5.8. Let X ~ SDS( 7 , A, k) with 0 < 7 < 1. Then 

E|Xf < 00 , for any 0 < r < 27 , 

E|X|'' = 00 , for any r > 27 . 


Proof. The moments of non-integer ord er ElX^ for any 0 < r < 2 can be computed using 
the following formula (see for example (IKlebanovl I2003L Lemma 2.2)): 


pOO 

E|xr = c, / (l-Re(/(<))) 

Jo 


dt 

r+i ’ 


with 

r 

r(l — r) cos(7rr/2) 

and where f{t) is the characteristic function of the distribution of X. Since SDS is 
a symmetric distribution, the characteristic function of X is real, and equal to 


f{t) = exp 



/ (1 - cos(f))(l + k) 
I — 2KCOs(t) -I- 1 



We may thus compute the moments. 


E|xr 



1 — exp 



/ (1 - cos(f))(l -I- k) 
i — 2Kcos(t) -1-1 



dt 


— Cf 



1 — exp 



/ (l-cos(f))(l-l-K) y] 
y — 2KCOs(t) -I- 1 j J 


dt 

^ 1+7 


~\~ Cr 



1 — exp 



y 1 - cos(t))(l -I- k) 

I _ 2Kcos(t) -I-1 



dt 

^T+7- 


Using the limit comparison test we see that the first integral converges for r < 27 and 
diverges for r > 27 , and the second integral converges for all r > 0 . □ 


5.5 Asymptotic behaviour 

In this Subsection we show that the tails of symmetric discrete stable SDS( 7 , A, k) distri¬ 
bution are indeed heavy with tail index 27 . 
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Proposition 5.9. The symmetric discrete stable distribution SDS( 7 , A, k) belongs to the 
domain of normal attraction of symmetric a-stable distribution with characteristic func¬ 
tion 


9W = exp{-|^AL±^li|2,|, 


Proof. Let Xi,X 2 , ■ ■ ■, Xn be i.i.d. SDS( 7 , A, k) random variables with characteristic func¬ 
tion 

(l — cos(t))(l + 


f{t) = exp < -A 


— 2 KCOs(t) - 1-1 


Let us denote Sn the normalized sum 

Xi X 2 ■ 


Sn = 


X„ 


7 , 1/27 


Then the characteristic function of Sn is given as 


E [«“•] = /" 


= exp < —A 


(1 — cos(t/ni/^'>'))(l -I- k ) 

k2 — 2«;cos(t/ni/27) -|- 1 


We use the Taylor expansion of cos to obtain 


logE le 


itSn 


\2ni/7 (1 — k)^ 


0{n 


- 3 / 27 )^ 


Hence 


=+ 3 / 27 )) 7 ^ as 77 -^ 00 . 

g{t) = lim E = exp i — ■;^ 7 r——. 

’ n-s-oo L i 27 (1 -k)27' ' j 


□ 


Theorem 5.10. Let X ^ SDS( 7 , A, k) with 0 < 7 < 1. Then 

- ^ 0 ^-^- 7/7^ i 


(5.3) 


hm cc^T-PdATl > x) = <{ r(i- 27 )cos(^ 7 )’ y ^ 2 ’ 

i 27 (1-k)27 t I */ 7 = 2 ■ 


Proof. We apply (|Ibragimov and Linnild . ll97lL Theorem 2.6.7.): SDS( 7 , A, k) distribution 
belongs to the domain of normal attraction of 8 ( 0 , /3,c,^) with a = 27 , fd = 0, c = 
A/2'>'(1 -I- k)'1'(1 — k )“27 and /i = 0, hence the tail functions of SDS( 7 , A, k) are given as 

F{x) = (ci-I-ai(x))|a:|““, for a: < 0, 

1 — F(a:) = (c 2 -I-a 2 (a^))a;““, for a; > 0, 

where ai{x) —>■ 0 as \x\ —>■ 00 . The constants ci,C 2 satisfy following conditions: 

/3 = (Ci - C 2 )/(Ci -I- C 2 ), 

r(l — a;)(ci-I-C 2 ) cos(7rQ:/2), if 
f (ci -I- C 2 ), 

We can easily see that for a y/ 1 we have 

1 A (l-hK)^ 


if 0 = 1. 


Cl = C2 = 


2 27 (1 - /c) 27 r(l - 27) cos(7r7) ^ 
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and for a = 1 we have 

A (1 + k)t' 1 

Hence 

lim x^'^f‘{\X\ > x) = lim x^'^{F{—x) + 1 — F{x)) 

x—^oo x—^oo 

= lim x^'^ [(ci + ai{-x))x~‘^^ + (c 2 + a 2 (a:))a;“^'’'l 

= 2ci. 

□ 


5.6 Asymptotic expansion of probabilities 

In this Subsection we give an asymptotic expansion of the probabilities of the symmet¬ 
ric discrete stable di stribution with k = 0. The f ollowing result is an adaptation of 
the approach used in Christoph and Schreibei ( 19981) for positive discrete stable random 
variables. 

Theorem 5.11. Let X ~ SDS( 7 , A), with 0 < 7 < 1. Then for any fixed integer m and 
n — 7 > 00 


m y-l-l 

(5.4) P(X = n) = --Tj- siiifjjTr)Bfyj + l,n — "fj) + 

77 - « ^ ^ ' 


f=i 


j! 


where B{x,y) = r(a:)r(?/)/r(a;-I-y) is the Beta function. Moreover 
(5.5) P(X = n) = - --- A'^r( 7 j -I- 1) 

'TT ‘ ^ 'I I 


t=i 


J'- 


Proof. Using the stochastic representation of S DS( 7 . A) random variable as a compound 
Poisson random variable with random intensity ( Slamova and Klebanov ( 2012h l we have 


POO 

P(X = n)= e~Mn{s)p^^{s)ds, 
Jo 


where In{s) is the modified Bessel function of the first kind and p^{s) is the density 
function of the random variable with characteristic function 

g{t) = exp{-A|t|'>'exp(-isgn(t) 77 r/ 2 )}. 


The density function pij{s) has the following series representation (jChristoph and Woli 

(Il992lb : 


(5.6) 


1 ( _1 \j-\-l 

p^is) = - V-n-A'’r(7j -h 1) sin(7j7r)s"'^-^"^ -k A^is), 

' TT ?! 

? = 1 


for any m > 0, where Amis) = 0{s 7 (m-i-i) i) ^3 g qo. We may compute the proba¬ 
bility as 


1 t-lV + l roo f-oo 

P(X = n) = - y^---A'^r( 7 j -h 1) sin( 7 j 7 r) / e“'*/„(s)s“'^''“Ms -k / e“®/„(s)Hm(s)ds. 

j! Jo Jo 


TT ^ -' 7! 

7 = 1 


38 

























We approximate the modified Bessel function /„(s) by the first term of its infinite series 
representation r(n + l)“^(s/2)". Then the first integral turns into 

The remainder term is obtained by computing the integral with j = m + I and by ap¬ 
proximating the ratio of two Gamma functions for large n using the Stirling’s formula 

(5.7) = n~'^^ (n~^ + O {n~^)) , as n —>■ oo. 

If we set m= [(7 + l)/7] and apply (15.71) on all terms in (15.41) . we obtain (15.51) . 

□ 

6 Properties of positive discrete stable random vari¬ 
ables with thinning operator of Chebyshev type 

The G thinning operator (of geometric type) used to define discrete stable distributions 
in the previous Sections is not the only possibility. As was showed in Chapter [5] we can 
consider also a T thinning operator (of Chebyshev type) given by the following probability 
generating function 


( 6 . 1 ) 


Q(z) = 




Ijm 


(1 + ^) + 

where p G (0,1), b G (—1,1) and m G N, and Tp{x) = cos (parccosx) . 

Theorem 6.1. A non-negative integer-valued random variable X is positive discrete stable 

with T thinning operator if and only if its probability generating function is given as 

( 6 . 2 ) 

{l + b)z^ 


V{z) = exp s — A arccos 


2 - (1 b)z^‘ 


with 7 G (0, 2], A > 0, 6 G (—1,1), m G N. 


Proof. Let h{z) = log'P{z). From Proposition 12.21 it follows that X is positive discrete 
stable if and only if h{z) = nh{Q{z)) for all n, where Q is as in (16.11) . Set 

,, , (l + 6)z--2&V 

h(z) = —A arccos —---- 

^ ^ V 2- (l-h6)z'" J 

and select 7 such that 1/p'’' = n. Then 

uni i f {l + b)Q{zr-2b^^ 

nh{Q{z)) = —An arccos 


= —An ( . 




arccos T, 


2-{l + b)Q{zy 
{l-^b)z^-2b 


= — An j p arccos 

= Hz)- 


^ ' 2-{l + b)z 
{l + b)z^-2b^'^ 
2 - (1 -h b)z^ 


□ 


We will denote the discrete stable distribution with Chebyshev thinning operator T 
and with parameters 7 G (0,2], A > 0, 6 G (—1,1) and m G N, by 7 ~PDS( 7 , A, b, m). If m 
is omitted then to = 1. If moreover b is omitted we will understand that 6 = 0. 
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6.1 Characterizations 

Theorem 6.2. Let 7 ' G (0,2] and 7 G (0,1] and assume that 7 ' < 27 . Let be a ^-stable 
random variable with Laplace transform exp{— m'*'}. Then 


rPDS( 7 ', A, b) = TPDS ( 777 , fc) . 

Proof. For sake of simplicity we will do the proof only for the case 6 = 0. The case 
670 can be proved in the same way. The probability generating function of X ~ 
TPDS ( 7 V 7 ) A^/'^'S.y) is computed as 


V{z) = F,z^ = Eexp J -A^/^S.^ ( . 


arccos■ 


2 - 


777 ' 


and using the Laplace transform formula for S.y we have 


V{z) = exp 



This is the probability generating function of TPDS( 7 ', A). 


□ 


Corollary 6.3. Let y,yi,l 2 , ■ ■ ■ be a sequence of i.i.d. random variables with probability 
generating function 


7?(z) = 1 — A arccos 


(1 + b)z — 2b 
2-{l + b)z' 


Let N be a random variable, independent of the sequence Yi,Y 2 ,.. ■, with Poisson dis¬ 
tribution with random intensity A^/'’' 7 rS.y, where 7 G (0,1] and S..y is a j-stable random 
variable with Laplace transform exp{— m'^}. A random variable X is positive discrete stable 
TPDS(7,A, 5) if and only if 

N 


Proof. Let X = Yj. Then T is a compound Poisson random variable with random 

intensity A^/'*'7rS.y and jumps • with characteristic function 


g{t) = 1 - ^ arccos 


(l + 5)e**-25 
2- (l + 5)ei* ■ 


The characteristic function of a compound Poisson random variable with intensity r 
and characteristic function of jumps h{t) is exp{— t(1 — h{t))}. Therefore X is in fact 
TPDS(1, 6)- We thus obtain the result from the previous Theorem 16.21 with 

7 ' = 7 . □ 


6.2 Continuous analogies 

Let us consider a T positive discrete stable random variable X ^ TPDS( 7 , A, 6 )- We are 
interested in the limit distribution of a random variable = aX, where a j, 0. We show 
that the limit distribution is in fact a-stable with index of stability 0 = 7/2 and with 
skewness 13=1. 
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Theorem 6.4. Let X be a random variable with probability generating function 

P( 2 ) = exp|-A |, 7 G (0,2], A > 0, 6 G (-1,1). 

Let X"' = aX and assume that X = Then the characteristic function of X°‘ converges 

pointwise to the characteristic function of a-stable distribution, 


rW=e,J-Afarcc„<^ + »“-^f 

^ ^ \ V 2 - (1 + 5)e“* 


J 07 TT'y f I+ b 

exp < —ct 2 ' cos — ' 


4 Vl-6 


7/2 


|t|T'/ 2 (l-isign(t) tan^ 


J 


Proof. For sake of simplicity we will do the proof only for 6 = 0. The characteristic 
function of can be approximated as 


oiat \ 'y 


2 - 6 “* 

1 + iatA 


log/“(t) = — A l^arccos 

« — A 1 arccos ■ 

\ 1 — iat 

Moreover arccos( 2 :) « — z as z —>■ 1. We have 

1 + iat —2iat 
1 — iat 1 — iat 


as a —>■ 0 . 


1 - 


Put together we obtain 


iogr(t)«- 


1.7/2 


—iat 
1 — iat 


as a —>■ 0 


—>■ —ct 2 ^(—it)'*'/^, as a —>■ 0 . 

Moreover we have (—it)^/^ = cos (l — isign(t) tan ^). The proof is therefore 

completed. □ 
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